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The dependence of the Josephson Plasma Resonance (JPR) frequency in heavily underdoped 
Bi2Sr2CaCu2 08+« on temperature and controlled pointlike disorder, introduced by high-energy 
electron irradiation, is cross-correlated and compared to the behavior of the afe-plane penetration 
depth. It is found that the zero temperature plasma frequency, representative of the superfluid 
component of the c-axis spectral weight, decreases proportionally with T c when the disorder is 
increased. The temperature dependence of the JPR frequency is the same for all disorder levels, 
including pristine crystals. The reduction of the c-axis superfluid density as function of disorder is 
accounted for by pair-breaking induced by impurity scattering in the CuC>2 planes, rather than by 
quantum fluctuations of the superconducting phase. The reduction of the c-axis superfluid density 
as function of temperature follows a T 2 -law and is accounted for by quasi-particle hopping through 
impurity induced interlayer states. 
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I. INTRODUCTION 



Significant controversy remains concerning an appro- 
priate model description of high temperature supercon- 
ducting cuprates (HTSC) in the underdoped regime, i.e. 
the regime in which the number of additional holes per 
Cu, p is smaller than the value 0.16 at which the critical 
temperature T c is maximum. 1 Whereas it is well estab- 
lished that the charge dynamics and transport proper- 
ties in the normal- and superconducting states in the 
overdoped regime (p > 0.16) are, by and large, deter- 
mined by well-defined quasi-particles, the role of quasi- 
particles in the underdoped regime is debated. The un- 
derdoped region of the (p, T) phase diagram is character- 
ized by several salient features^ At T* > T c , the well- 
known "pseudo-gap" in the excitation spectrum opens 
up. This has been interpreted as either being related 
to the advent of another type of (spin- or charge-) or- 
der competing with superconductivity, driving T c down 
as p is diminished, or, alternatively, as signalling the 
formation of precursor Cooper pairs without long-range 
phase coherence. Then, T c is interpreted as the demise 
of long range superconducting phase order due to strong 
thermal 3 -"^ or quantum 6 -^ phase fluctuations. Strong 



support for this scenario has come from the violation 
of the Glover-Tinkham-Ferrell conductivity sum-rule ap- 
plied to the c-axis spectral weight^ also, the linear re- 
lation between T c and the superfluid density 8 has been 
interpreted as the result of T c being determined by phase 
fluctuations in a Kosterlitz-Thouless-type scenario i^Jil A 
smoking gun for such a scenario would be an important 
reduction of the c-axis superfluid density p c s , which is 
determined by Cooper pair and quasiparticle tunneling 
between adjacent strongly superonducting Cu02 layers 
through the weakly superconducting rocksalt-like block- 
ing layers, with respect to the in-plane stiffness in 
underdoped cuprates below T c . 

However, apart from phase fluctuations, other mech- 
anisms for the reduction of p c s , arising from disorder in 
the crystalline structure of underdoped cuprates, can- 
not be ignored*^ First, scanning tunneling spectroscopy 
(STS) experiment o 12 ! 13 ' 14 have revealed large variations 
of the magnitude of the gap maximum A pea k, as inter- 
preted from conductance curves measured on the surface 
of P^S^CaC^Os+s crystals. This has motivated recent 
interpretations of weakened c-axis superfluid response in 
this materia l 15 ' 16 as well as in La2- 3; Sr 3 ;Cu04 17 in terms 
of finely dispersed 5-20 nm-sized non-superconducting 
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regions within the C11O2 planes. Such regions may arise 
from the suppression of the superconducting order pa- 
rameter by dopant atoms p» such as out-of-plane oxygen 
atoms in the Bi2Sr2CaCu2 0s+5 compound^ 

Moreover, the d-wave symmetry of the gap function is 
at the basis of several mechanisms by which pointlike dis- 
order reduces the c-axis superfluid density. The appear- 
ance of quasiparticle (virtual-) bound states and their 
smearing by a finite defect density leads to an increase of 
the density of states (DOS) near the nodal (it, ir) direc- 
tions (the so-called lifetime effect) The same pointlike 
defects increase the quasiparticle scattering rate T s . It 
was conjectured that in the case of coherent (in-plane 
momentum preserving) quasiparticle tunneling, the can- 
cellation of these effects leads to disorder-independent 
low-temperature c-axis quasiparticle conductivity and 
superfluid density^ The approach of Ref. [2l[ neglects 
the crystal structure of the tetragonal HTSC,— which 
leads to the dependence of the interlayer hopping inte- 
gral t± on the in-plane momentum (k x ,k y ). For sim- 
ple tetragonal structures, interlayer hopping occurs via 
Cu 4s orbitals in adjacent planes. Its momentum de- 
pendence t± = ^[cosfc^a — cosfcj,a] 2 is determined by 
the in-plane overlap of the bonding oxygen 2p level with 
the 4s level of the neighboring Cu atom i 22 i 23 As a re- 
sult, c-axis tunneling occurs nearly exclusively for the 
anti-nodal directions at which quasiparticles are unlikely 
to be excited. In body-centered tetragonal structures 
such as Bi2Sr2CaCu208+5, hopping is also suppressed 
along the (k x ,k y ) = (tt,0) and (0,7r) lines, yielding 
t± = t° ± [cos k x a — cos k y a] 2 cos ^fc^acos \k v a^ In either 
case, disorder is always relevant for the nodal directions. 
Then, from the lifetime effect, one expects a quadratic 
decrease with temperature of the reduced c-axis super- 
fluid density-2&2£ 



strong forward scattering, the result 
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Here a c is a dimensionless constant of order unity and 
the parameter A was assumed, in Refs. [23j and |25|, to 
correspond to the maximum amplitude of the Bardeen- 
Cooper-Schrieffer d-wave gap. Eq. (fl]) essentially differs 
from that derived for the a6-plane superfluid density 
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in that the leading temperature-dependent term has a co- 
efficient that is smaller by a factor (r s /Ao) 2 ^ The pres- 
ence of defects in the rocksalt-like (BiO) layers tends to 
break the d-wave symmetry of the hopping integral, and 
renders quasiparticle hopping possible for other values of 
the in-plane momentum, and notably along the order pa- 
rameter nodes i 11 ' 23 ' 25 ' 26 A condition for this "impurity- 
assisted hopping" (IAH) to be effective is an anisotropic 
scattering matrix of the interplane defects. Notably, for 
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was obtained 
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for T s < k B T < 
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T C 2&ZL Here V x is the 
magnitude of the impurity scattering potential of the 
out-of-plane defects and N(Ep) is the density of states 
at the Fermi level in the normal state. The effect of im- 
purities can be distinguished from that of boson-assisted 
interlayer hopping; for the latter, a very similar result is 
obtained, but with the leading temperature-dependent 
term proportional to T 3 . 25 ' 28 Finally, direct hopping of 
quasiparticles was suggested to lead to a small, linearly 
temperature-dependent, reduction of p c s *^ 

In this paper, we address the mechanism by 
which the c-axis superfluid density in underdoped 
Bi2Sr2CaCu208+,5 ( with p = 0.1Q ) is reduced by us- 
ing disorder, in the form of Frenkel pairs introduced by 
high energy electron irradiation, as an independent con- 
trol parameter. Electron irradiation, the effects of which 
are taken to be similar to those of Zn-doping, 2 — has pre- 
viously been used to study the effect of pointlike disorder 
on the resistivity, critical temperature? 3 ^ and Nernst ef- 
fect of YBa 2 Cu 3 07 and YBa 2 Cu 3 06.6'^ In the latter 
material, electron irradiation eventually leads (at high 
fluences) to the breakdown of the well-known Abrikosov 
Gor'kov relation 3 ^ 3 - 
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(with T c0 the critical temperature when the normal state 
scattering rate T is equal to zero, and ^ the digamma 
function) as well as a significant increase of the fluctua- 
tion regime near T c £^ Both effects were interpreted as the 
effect of strong superconducting phase fluctuations . 30 ! 31 
The in-plane and c-axis superfluid densities of Zn-doped 
YBa2Cu307-«5 were studied by Panagopoulos et alM- and 
by Fukuzumi, Mizuhashi, and Uchida.— The progressive 
inclusion of Zn leads to a rapid decrease of the in-plane 
superfluid density corresponding to an increase of 
the in-plane penetration depth A a {,(0) oc (p° i b )~ 1 t' 2 , and 
a more modest decrease of p% oc cr c (T c ) oc T c , that vio- 
lates the c-axis conductivity sum-rule^^ [a c {T c ) is the 
c-axis conductivity at T c ]. As for the low-T tempera- 
ture dependence, a gradual change of both p^ b and p c s 
from T-linear to T-squared has been reported.— Stud- 
ies on Bi2Sr2CaCu208+<5 are limited to electron irra- 
diation of the single crystalline optimally doped ma- 
terial, that show a linear decrease of T c with electron 
fluencei 36 ' 37 ! 38 The c-axis superfluid density in a under- 
doped pristine Bi 2 Sr2CaCu208+i single crystal has been 
previously studied by Gaifullin et ai, who invoked the 
IAH model to explain the much stronger temperature de- 
pendence of p c s in underdoped with respect to optimally 
doped Bi 2 Sr 2 CaCu208+5^ 
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Below, we report on c-axis coupling in the supercon- 
ducting state measured through the Josephson Plasma 
Resonance (JPR) ) 40 ' 41 ' 42 ' 43 which, in our underdoped 
Bi2Sr2CaCu208+5 crystals takes place in the microwave 
frequency regime below 70 GHz. The JPR frequency 
f p i is sensitive to the value of Ao, as well as to fluctu- 
ations of the superconducting order parameter phase in 
the Cu02 planes^ The evolution of f p i(T) with temper- 
ature depends simultaneously on the quasiparticle dy- 
namics and on the strength of fluctuations; the plasma 
resonance peak is broadened both by the quasiparticle 
tunneling rate and by crystalline disorder. 17 ' 44 However, 
the dependence of / 2 ; oc p c s on the disorder strength is ex- 
pected to be quite different, depending on which mech- 
anism is predominant. In the following, we show that 
the disorder dependence of the c-axis plasma frequency 
is a sensitive probe, that allows one to identify in detail 
what physical mechanisms are at the basis of the reduc- 
tion of the superfluid density in Bi 2 Sr 2 CaCu2 08+5. It 
turns out that, even in our heavily underdoped crystals, 
(incoherent) c-axis quasiparticle hopping is essential for 
a consistent description of the data. We find that the en- 
ergy scale Ao, which turns out to be Deltag ~ 2.5fcsT c 
for all underdoped crystals, is to be interpreted as an 
energy scale governing nodal quasi-particle excitations. 




T(K) 

FIG. 2: (color online) Real ( R s ) and imaginary ( X s ) parts 
of the surface impedance Z s of underdoped Bi2Sr2CaCu20s+5 
crystals, irradiated with 7.7x 10 18 and 6x 10 19 electrons cm -2 , 
respectively. The data were obtained from the resonance fre- 
quency shift and the quality factor of a superconducting Pb 
cavity operated in the TEon mode. The inset shows the ab- 
plane- and the c-axis dc resistivity of the crystal irradiated 
with 6 x 10 19 electrons cm -2 . 



II. EXPERIMENTAL DETAILS 

The underdoped (T c = 65 ± 0.5 K, p » 0.10) 
Bi 2 Sr2CaCu 2 08+<5 single crystals, of typical dimensions 
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FIG. 1: (color online) Real part of the ac magnetic sus- 
ceptibility of the investigated underdoped Bi2Sr2CaCu2 08+a 
crystals. The crystals were irradiated at 22 K with 2.3 MeV 
electrons to the indicated fluences. The ac field amplitude 
h ac = 4.2 mOe, the ac frequency was 560 Hz. The curves 
show a screening onset determined by a surface layer con- 
taining optimally doped material; this screening is suppressed 
when hac > 0.5 Oe. The steep drop at the lower end of the 
transitions corresponds to bulk screening by the underdoped 
material. 



500 x 300 x 40/mi 3 , were selected from the same boule, 
grown by the travelling solvent floating zone method at 
the FOM-ALMOS center, the Netherlands, in 25 mBar 
O2 partial pressure. 4 - The crystals were annealed for one 
week in flowing N2 gas. We have also measured a set of 
optimally doped control samples ( T c = 86 K ). These 
were also grown by the travelling solvent floating zone 
technique, at 200 mbar oxygen partial pressure, and sub- 
sequently annealed in air at 800°C. The crystals were ir- 
radiated with 2.3 MeV electrons using the Van de Graaff 
accelerator at the Laboratoire des Solides Irradies. The 
beam was directed along the crystalline c-axis during the 
irradiation. To prevent recombination and clustering of 
point defects, the irradiation is carried out with the crys- 
tals immersed in a liquid hydrogen bath ( 22 K ). The 
electron flux is limited to 2 x 10 14 e"cm" 2 per second. 
Crystals UD5-UD8 were irradiated to a total fluence of 
0.53 x 10 18 , 3 x 10 18 , 7.7 x 10 18 , and 8.8 x 10 19 e'cnr" 2 re- 
spectively. After measurements, crystal UD5 was irradi- 
ated a second time to a total fluence of 6.0 x 10 19 e~cm -2 
and was henceforth labeled UD5b. The high energy elec- 
tron irradiation creates random atomic displacements in 
the form of Frcnkcl pairs, both in the Cu02 bilayers and 
in the intermediate cation layers, throughout the sam- 
ples. 

The superconducting transition temperature T c was 
determined by ac susceptibility measurements using a 
driving field of amplitude 4.2 mOe and a frequency of 
560 Hz, directed parallel to the c-axis. For all under- 
doped crystals, the superconducting transition is rather 
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FIG. 3: Microwave absorption, measured using the TMnio 
mode ( 19.2 GHz ) of one of the OFHC copper cavities, of 
three of the underdoped Bi2Sr2CaCu20g+5 crystals before 
and after irradiation with 5.3 x 10 17 , 3 x 10 18 , and 7.7 x 10 18 
2.3 MeV electrons cm -2 , respectively. 



broad. The superconducting transition takes place in two 
steps: there is a slow increase of screening at high temper- 
ature, followed by a rapid step of the diamagnetic signal 
at lower temperature ( Figure [T]). The high temperature 
screening vanishes when the excitation field amplitude is 
increased beyond 0.5 Oe, while the step at lower tempera- 
ture is robust. This shows that doping is macroscopically 
inhomogeneous, and that the crystals are surrounded by 
a thin surface layer of higher doping. This layer could 
not be eliminated by cutting the crystals. The overall 
shape of the transition is unaffected by the electron ir- 
radiation. The transition widths are of the order of 4 
K, which is usual for such low doping. After irradiation, 



the transition widths slightly increase. For all crystals, 
the transition to zero dc resistivity and bulk supercon- 
ductivity occurs at the temperature at which the lower 
screening step takes place, see e.g. the Inset to Fig. [2] 
Therefore, the lower temperature feature was adopted as 
characterizing the bulk T c of the underdoped crystals. 

Crystals were further characterized by the measure- 
ment of the temperature variation of the in-plane pene- 
tration depth, A a6 (T)/A ah (0) - 1 = AA a6 /A a6 . For this, 
a crystal is mounted on a sapphire rod, in the center of a 
superconducting (Pb) resonant cavity immersed in liquid 
4 He, and operated in the TE n mode. The cavity reso- 
nant frequency was / ~ 27.8 GHz, and the quality factor 
Q ~ 4 x 10 5 . The crystal is mounted in such a way that 
the magnetic microwave field is perpendicular to its ab 
plane and solely in-plane screening currents arc induced. 
From the shift A/ of the cavity resonance frequency in- 
duced by the sample, we determine the surface reactance 
X s = 2n^ G 2 Af (with p = 4tt x 1(T 7 Hm" 1 ). The sur- 
face resistance R s — 2ir[ioGi(AQ)~ 1 was obtained from 
the change of the quality factor. The geometrical fac- 
tors G\ and Gi were determined by comparing the sur- 
face impedance in the normal state, X s = R s = irfxof^sj 
to the value expected from the normal state resistivity, 
p = 7r hq f '81 . 46 It was retrospectively checked that all 
measurements were carried out in the skin effect regime, 
in which the normal state skin depth 5 S is much smaller 
than the sample dimensions. The relative change of 
the penetration depth was determined from the behav- 
ior of the surface reactance at low temperature, at which 

X s P3 27T/i /A a 6. 

Figure [2] shows that the temperature at which the de- 
crease of the surface resistance R s was observed corre- 
sponds to the (high temperature) onset of screening in 
the ac susceptibility measurement. This indicates that 
the surface skin depth of thickness ~ 7/^m probed by the 
microwave field contains patches with larger hole content 
P- 

The JPR measurements were performed using the cav- 
ity perturbation technique, using the TMoi n modes (n — 
0, ...,4) of two Oxygen-Free High Conductivity (OFHC) 
copper resonant cavities ( Q ~ 4000 — 10000 ), mounted 
on a cryocooler cold head. The measurement frequencies 
ranged between 19.2 and 39.6 GHz.— Further measure- 
ments were made applying the bolometric technique, us- 
ing waveguides in the TEoi travelling wave mode^ In 
both measurement set-ups, the electrical microwave field 
is applied along the c-axis of the crystal. In contrast 
to the previously described surface impedance measure- 
ments, screening of the electric microwave field is very 
poor because of the high electronic anisotropy of the 
crystals. The underdoped samples are in the complete 
depolarization regime and thus the bulk electromagnetic 
response is probed. By monitoring the power absorption 
as a function or frequency for a fixed temperature, the 
Josephson Plasma Resonance is detected as a sharp ab- 
sorption peak in the microwave response, see Fig. [3] We 
determine the JPR frequency at a given temperature, 
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f p l(T), as the measurement frequency at the tempera- 
ture at which dissipation is maximum. 



III. RESULTS 

Figure [4] collects the values of the critical temperature 
as function of electron dose, for the set of underdoped 
samples as well as the optimally doped control samples. 
Both the T c of the underdoped and the optimally doped 
crystals decrease linearly with irradiation fluence. The 
derivative of T c with respect to fluence of the optimally 
doped crystals concurs with that measured by Behnia et 
al£l and Nakamae et al£2- but is two times lower than 
that measured by Rullier-Albenque et al..— The overlap 
between the variation with fluence of the screening onset 
temperature in the underdoped crystals and the T c of the 
optimally doped crystals shows that the thin surface layer 
on the underdoped Bi 2 Sr2CaCu 2 08 + 5 hasp ~ 0.16. The 
critical temperatures, normalized to the critical temper- 
atures T c q of the unirradiated crystals, can be superim- 
posed on Eq. ((4]), yielding estimates of the normal state 
scattering rate T (see Fig. [5]) . This procedure supposes 
that T c q corresponds to the critical temperature in the 
absence of disorder; we shall see below that this is not 
justified, so that the estimated T values are in fact lower 
limits for each crystal. 

The relative change with temperature of the in-plane 
penetration depth X a b is depicted in Fig. [6j For all 
underdoped crystals, including the unirradiated ones, 
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FIG. 4: Variation of the critical temperature with electron 
fluence for single-crystalline underdoped Bi2Sr2CaCu2 0s+,5 
(•). A comparison with the T c of a series of optimally doped 
control crystals (□) shows that the onset temperature of mag- 
netic screening (o) corresponds to flux exclusion by a surface 
region containing optimally doped material. 
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FIG. 5: Critical temperatures T c normalized to the critical 
temperature T c o of the unirradiated crystals, and superim- 
posed on the Abrikosov-Gor'kov relation @. Note that the 
actual data points should be shifted downwards along the 
curve, because T c o does not truly correspond to the critical 
temperature in the absence of disorder. 




500 1000 1500 2000 2500 3000 3500 

T 2 ( K 2 ) 

FIG. 6: (color online) Relative change of the afc-plane pene- 
tration depth A\ ab /\ ab (0) = [X S (T)/X S (T -> 0)- l]/27r/i /, 
as obtained from the change in the inductive part X 3 of the 
surface impedance. 



^ab(T) — A a h(0) varies quadratically with temperature 
at low T. Such a temperature dependence has been 
associated with quasiparticle scattering in the unitary 
limit by point defects situated within the CuC>2 planes 
of the d-w&ve superconductor, i.e. N(Ef)V 3> 1, T ~ 
n d /nN(E F ), and T s ~ 0.6(rA ) 1/2 4£ Here V is the scat- 
tering potential of the defects in the planes, with density 
n<j. The quadratic temperature dependence of AA Q {,(T) 
is at odds with a possible important role of thermal 
phase fluctuations, for which a linear T-dependence was 
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FIG. 7: Temperature dependence of the JPR frequency of 
electron-irradiated underdoped Bi2Sr2CaCu20g+ $. 



predicted.— As for the magnitude of the T 2 -contribution 
to A a h, a very modest change is found for the lower irradi- 
ation fiuences. Only for fluences exceeding 10 19 e _ cm~ 2 
does the in-plane penetration depth increase significantly 
with defect density. 

We now switch to the central results concerning the 
Josephson Plasma Resonance. Figure [7] shows the JPR 
frequency f v \ (T) of crystals UD5-UD7 as function of tem- 
perature ( measured in Earth's magnetic field ). The 
temperature at which f p i (T) extrapolates to zero is well- 
defined and corresponds to the critical temperature of 
the bulk, underdoped portion of the crystals, i.e. the 
main transition in the ac susceptibility and zero resis- 
tance. This shows that the JPR probes the c-axis re- 
sponse in the heavily underdoped bulk and is insensi- 
tive to the surface quality of the samples. From Fig. 
one sees that not only T c , but also f p i is strongly de- 
pressed by the electron irradiation. Fig. [8] collects val- 
ues of the low-temperature extrapolated value f p i (0) ver- 
sus the critical temperature, and reveals the propor- 
tionality between / 2 ; (0) and T c . This dependence is 
clearly different from the variation of f p i with oxygen 
doping. The same Figure recapitulates results for dop- 
ing levels p = 0.1322- and O.ll^ 7 - the evolution of which 
recalls the exponential /p,(0)(T c )-dependence found by 
Shibauchi et al.. le The results are somewhat similar to 
those obtained by Fukuzumi et al. for Zn-doped and 
oxygen-deficient YBa2Cu307_a: 3 » however, their Fig. 5 
shows a linear dependence of f p i on T c for underdoped 
YBa2(Cui_ a; Zn a ;)306.63, and a hyperbolic or exponential 
evolution of f p i{T c ) as function of 5. 

As for the temperature variation of f pl (T) oc p c s , it 
turns out to be identical for all the underdoped crystals, 
including the pristine ones, and does not depend on the 
defect density. Fig. © shows f p i(T)/f^{0) plotted ver- 
sus the reduced temperature T/T c . For hf p i < T, this is 
representative of the c-axis superfluid fraction p c s . The 
same graph may also be interpreted as the maximum 
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FIG. 8: /pi(0) versus T c (open symbols), as well as /„i(0) oc 
Ps(0) oc jc(0) versus T c (closed symbols), for both the under- 
doped irradiated crystals (squares) and a set of crystals with 



different doping levels (p s 



i 0.13,— and p ■ 



0.11 4 '). 



Josephson (c-axis) critical current, 
f c = 2^h/2e)ee s- 1 f pl = 



AirepaX^s 



(5) 



normalized to its value for T — * 0. Here s is the spacing 
between Cu02 planes, the c-axis dielectric constant e ~ 
11.5, 5 & and e = 8.854 x 1CT 12 Fm" 1 . 



IV. DISCUSSION 

To distinguish between the different mechanisms re- 
sponsible for the reduction of the c-axis plasma fre- 
quency as temperature is increased, we dispose of three 
tools. First, there is the variation of f p i(0) with disorder 
strength, which manifests itself starting from the smallest 
electron fluences, and proportionally follows the evolu- 
tion of T c with disorder. Second is the temperature vari- 
ation /pi(T)// 2 z (0), that follows a 1 - a(T/T c ) a depen- 
dence, with a i~ 2 independent of the disorder strength. 
Finally, there is the comparison with the behavior of the 
in-plane penetration depth, X a b(T) / X ab (0) ~ 1 + f3T 2 . A 
successful model description should account for all three 
dependences correctly. 

The theoretically expected low-temperature depen- 
dence A C (T) depends strongly on a number of circum- 
stances. First is the question whether superconductive 
coupling is three-dimensional ( i.e. the c-axis momen- 
tum k z is a good quantum number ), 23 or whether it 
is mediated by Josephson tunneling between CuC>2 lay- 
ers. Josephson coupling can be weakened by direc t 21 ' 27 or 
boson-assisted quasiparticle tunneling ) 25 ' 27 ' 28 or by tun- 
neling that involves intermediate defect-induced states 
between the layers (IAH) i 23 ' 27 For "direct" tunneling, 
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e.g. occurring through direct overlap of the supercon- 
ducting wave functions in the CuC>2 planes, one dis- 
tinguishes the case of conserved in-plane momentum k\\ 
( "coherent" tunneling) from the case where it is not pre- 
served ( "incoherent" tunneling 2 ^ — this situation would 
yield a vanishing j£ in a ci-wave superconductor). A pre- 
vailing effect of nodal quasiparticles leads to a more rapid 
decrease of j° , see Eq. (fTJ). Finally, in all cases the 
anisotropy of the transfer matrix t± is expected to have 
an important influence on the temperature dependence 
of A e ) 23 i 52 notably reconciling a weak T-dependence of A c 
with a strong variation of X a b(T). Here, the experimen- 
tally observed temperature dependence of /p;(T)// 2 ; (0) 
actually allows one to discard a dominant role of a pos- 
sible d-wave symmetry of the transfer matrix t± j22 since 
( for ksT ^> hT s ) this leads to a weak temperature de- 
pendence, fpi(T)/f^(0) ~ 1 - aT 5 , observed in slightly 
overdoped Bi 2 Sr 2 CaCu208+ j 39l5 ° and optimally doped 
HgBa2Ce n _iCu„02n+2+(Sr^ but not in the present data 
on underdoped Bi2Sr 2 CaCu 2 08+,5. The modification of 
ij_ arising from the body-centered Bravais lattice of the 
Bi 2 Sr 2 CaCu208+,5 compound will influence the maxi- 
mum c-axis critical current. However, it will not change 
the expected 1 — aT 5 temperature dependence, since this 
arises from the specific coincidence in fc-space of the zero 
of t±(k x ,k y ) with the nodal direction of the order pa- 
rameter. Thus, models for superconductive coupling, 23 
or direct Josephson coupling with a vanishing hopping 
integral along the nodal line^ are in inadequacy with 
the data. 

We next exclude a dominant role of direct quasiparti- 
cle tunneling. Even though the similar T 2 -dependences 
of the low-temperature a6-plane and c-axis penetration 
depths suggests such coupling, the disorder dependence 
is at odds with the experimental result. Radtke et al^ 
and Latyshev et al^ find that for direct coupling, the 
low temperature c-axis critical current 



3 



c.direct 



7ro£(0)Ao _ 4Tret 2 ± N n (E F ) 



(6) 



is, for ksT <C HT S ^ 20 — 30 K, to lowest order indepen- 
dent of the defect density due to the cancellation of the 
scattering-rate dependences of the quasi-particlc conduc- 
tivity a c q and Aq. The model was further worked out by 
Kim and Carbotte, who find that to first order 



•c, direct -i 



VTITA 2 



1 -a- 



1 



Vi + Ao/r 



(7) 



both for the case of constant t± (where a « 1) and 
angular-dependent tj_ (a = if )^ The nonlinear depen- 
dence ([7]) is at odds with the observed linear evolution of 
j£ with irradiation fluence. 

The temperature and disorder-dependence of the low 
temperature c-axis JPR frequency is more successfully 
described by a model for incoherent tunneling. According 
to Latyshev et ai, an incoherent tunneling process yields 
jc,i _ jc,direct Ao / Ep _ 4Tret 2 x N n (E F )A /hE F . The ex- 
tra factor Aq then explains the linear relation between 
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FIG. 9: Square of the JPR frequency, normalized to its 
low temperature extrapolation f P i{0), versus reduced tem- 
perature T/T c . This plot is representative of the tempera- 
ture dependence of the c-axis superfluid stiffness, as well as 
of the maximum c-axis Josephson current: fpi{T)/fpi(0) ~ 
A»(0)/A2(T) ~ p c s {T)/pl{Q) ~ f c (T)/f c {Q). The drawn line is 
a fit to Eq. © with A = 2.5k B T c . 



jc(0) and T c , accepting that in a d-wave superconductor 
with impurity scattering in the unitary limit, Ao(r) is 
simply proportional to T c (r). 33 The linear dependence 
on A is found in the 1AH model, see Eq. (gj) .23,27,28 
The latter expression consistently describes the fact that 
the temperature dependence of f p i does not change with 
defect density: the temperature dependent term writes 
(T/Ao) 2 oc (T/T c ) 2 . For the same reason, the "lifetime 
effect", Eq. dTJ, does not describe the data: in the uni- 
tary limit, the leading temperature-dependent term has 
an extra factor T s /Aq ~ (r/Ao) 1 / 2 and is therefore ex- 
pected to strongly depend on defect density. The ob- 
served defect-density independence of p c s (T)/p c g would 
require the strength of the scattering potential of the in- 
dividual irradiation defects in the Cu0 2 planes to be in 
the Born limit, which contradicts the results on the tem- 
perature dependence of the a6-plane penetration depth. 
We note that the toy model for incoherent hopping of 
Ref. [IH, which yields A" 2 oc [1 - A(r s /A ) 2 - . . .] ~ 



1 



also describes the initial linear de- 



&(r/Ao)-. 

crease of j£(T c ) (Eq. (31) of Ref. [ll|). 

Given the success of the IAH-model in qualitatively 
explaining the temperature- as well as the disorder de- 
pendence of the JPR data, we perform a direct fit of 
/ 2 Z (T)// 2 (0) to Eq. ©, shown in Fig. M The only 
parameter is the value Ao = 2.bksT c . This value not 
only means that T s > 30 K, comforting our interpre- 
tation of the decrease of j% oc T c in terms of unitary 
scatterers induced in the Cu02 planes, it is also remark- 
ably close to characteristic energy scales found in recent 
Raman scattering^ and STM experiments^ The first 
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study finds that in underdoped HgBa2Cu04, the Bi g 
Raman mode, which directly couples to the same low- 
energy nodal quasi-particle excitations that are responsi- 
ble for the reduction of the c-axis superfluid density in the 
IAH model, is characterized by an energy scale ~ 2ksT c 
(related to the nodal slope of the gap function). 53 We 
surmise that Ao is a closely related parameter. The sec- 
ond study finds that the total tunneling gap amplitude 
is determined by two energy scales, the smaller of which, 
w 2.8/cbT c , is related to superconductivity— 4 It is 
interesting that the existence of a second small energy 
scale describing nodal quasiparticle exitations may well 
be responsible for the observed suppression of j° as func- 
tion of doping (decreasing p or S) The proportionality 
of this decrease to the ratio of the gap maximum (at the 
antinode) and T c finds a logical explanation if the smaller 
energy scale ( ~ 2 — 2.5fcsT c ) is what determines the 
magnitude of the c-axis critical current density. 

Thus, the (incoherent) interlayer assisted hopping is 
a feasible candidate for the reduction of the c-axis su- 
perfluid density in underdoped Bi2Sr2CaCu20s+(5: it 
parametrically describes the data, and numerical val- 
ues extracted for the relevant energy scale determining 
the quasiparticle excitations is the same as found in 
other experiments. The model does have several caveats: 
first, there is the above-mentioned puzzle that it requires 
the temperature dependence of the afe-plane superfluid 
density to be explained by the lifetime effect, whereas 
the same effect does not seem to play a role in the c- 
axis superfluid density, other than providing the quasi- 
particles. The second is the identification of the defects 
in the rocksalt-like layers that are responsible for inter- 
layer scattering. The model by Xiang et a/i 23 ' 25 requires 
these out-of-plane defects to be weakly scattering, with 
a strongly anisotropic potential that leaves the reflected 
wavevector close to the incident ( "strong forward scatter- 
ing"). Although candidates may be out-of-plane oxygen 
defects 1 ^ or cation disorder, the constraints imposed on 
the scattering potential seem very strict. In the end, 
the agreement of defect-density independence of the ex- 
perimental fpi(T)/f p i(0) with the IAH prediction^ 3 - may 
be completely fortuitous. The coincidence of the tem- 
perature dependence of f p i of the pristine crystals with 
even the most heavily irradiated ones indicates that dis- 
order plays an important role in all samples. Notably, 
we expect T c of pristine crystals to be substantially sup- 
pressed with respect to that of hypothetically "clean" 
underdoped Bi2Sr2CaCu208+<s. Furthermore, impurity 
scattering is likely to completely suppress any role of 
quasi-particles in the c-axis electromagnetic response of 
this compound, leaving only pair tunnelingi&S 5 - 

The remaining mechanism for the reduction of the su- 
perfluid density in the presence of pair tunneling only is 
that of order parameter phase fluctuations in the CuC>2- 
layers. The effect of quantum phase fluctuations on the 
a6-plane and c-axis superfluid densities was examined by 
Paramekanti et al.fi^ who performed an analytical study 
of an XY model for the superconducting order phase <f> 



on a two-dimensional (2D) lattice of spacing £o (repre- 
senting the coherence length), within the self-consistent 
harmonic approximation. The authors conclude that, in 
contrast to thermal phase fluctuations, quantum phase 
fluctuations in quasi 2D high temperature superconduc- 
tors are important at all temperatures. The low carrier 
density in these materials leads to inefficient screening of 
the Coulomb interaction between charge carriers, and a 
sizeable reduction of the magnitude of the a5-plane su- 
perfluid density (without change of its temperature de- 
pendence). Within this model, the JPR frequency is also 
renormalized downwards because of fluctuations of the 
phase in the layers: 

fpi (T) = #(0)e-<**> » #(0) (l - |<0i) -...). 

(8) 

Paramekanti et al. estimate the phase difference between 
two points separated by a vector perpendicular to the 
superconducting layers as {(f>\_) w v /e 2 /47reeoCo£o(0)s, 7 
with £ S = h 2 s/16e 2 TTfioXl b (0) the in-plane phase stiff- 
ness. We observe that, given the Uemura relation s s oc 
fceT c ., 8 the resulting expression naturally describes the 
experimentally observed exponential depression of j° as 
function of doping. 1 ^. However, even if one explicitly de- 
velops the dependence of eo-s in terms of the variance of 
the in-plane phase (</)|) ~ {<j)\), Eq. © fails to describe 

the linear / 2 z (T c )-dependence (Fig. [8]). Therefore, the re- 
duction of p c s with increasing disorder cannot be ascribed 
to quantum phase fluctuations only - pair-breaking in the 
Cu02 layers must play a significant role. 

A noteworthy prediction of the quantum fluctuation 
scenario is that the temperature evolution of the c- 
axis superfluid density is entirely determined by that 
of the in-plane phase stiffness, i.e. the c-axis and ab- 
plane superfluid densities follow the same dependence 
at low temperature. Inserting the experimental result 
A Q 6 = A a6 (0)(l + /3T 2 ) into the prediction of Ref. 0], one 
would expect 

d[f P i(T)/f P i(0)] 2 = _Ci , T 2 / 2^ 

d{T/T c y 4 P c \/4^eo£oeo(0) S ' K ) 

This is experimentally verified; taking the data of Fig. [6] 
and A afc (0) « 300 nm 4 ! we find that Eq. © is obeyed 
with Ci « 0.3 ( C\ should be order unit}*£ ). The ex- 
perimental independence of f p i (T) / f p i (0) on defect den- 
sity demands that /3T 2 A a ;,(0) is disorder independent. 
Adopting the generally accepted view that A~ b 2 (T) is de- 
scribed by Eq. ^ with T s in the unitary limit, this would 
imply that /?T 2 A ab (0) ra a ah \ ab (0){T c /T)^ 2 (T c /A p 2 , 
and therefore that A Q b(0) -2 oc T c /T. This is as yet un- 
verified, as the different sizes and aspect ratios of our 
crystals prohibit a direct comparison of the absolute val- 
ues of X a b- 

Note that the case of screening by nodal quasiparticles 
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was also studied in Ref. Q. Then, 

d\f P i(T)/f pl (p)]\ WT 2 C 

d(T/T c y ~ a q ' 1 ' 

where W q is the a6-plane quasiparticle sheet conductiv- 
ity, normalized to the quantum conductivity e 2 /h. This 
formula also describes the temperature dependence of 
the data, provided that a q « 3; moreover, the ratio 
j3 l~u q ~ a ah A (T c /A o )m/iV s (0)e 2 should be disorder- 
independent (N S (Q) is the quasiparticle density of states 
and m the effective mass). Given the theoretical expec- 
tation N S (Q) ~ T 1 / 2 [48] and A oc 1 - T [33], this model 
again fails to describe the reduction of the zero temper- 
ature c-axis superfluid density in terms of quantum fluc- 
tuations only. 

V. SUMMARY AND CONCLUSIONS 

We have cross-correlated the dependence of the c-axis 
Josephson Plasma Resonance frequency in heavily under- 
doped Bi2Sr2CaCu20s+5 on temperature and controlled 
disorder (introduced by high energy electron irradiation) , 
and compared both with the behavior of the in-plane 
penetration depth. It is found that the c-axis critical 
current is depressed with increasing disorder strength, 
proportionally to the critical temperature T c . Both the 
in-plane- and out-of-plane superfluid densities follow a 
T 2 temperature dependence at low T. The tempera- 
ture dependence of the c-axis response is independent 
of disorder, indicating that we are probing the super- 
fluid density. The superfluid response of the pristine un- 
derdoped crystals is indistinguishable from that of heav- 
ily irradiated ones, suggesting that pristine underdoped 
Bi2Sr2CaCu208+,5 commonly contains sufficient disorder 
in the Cu02 planes for the critical temperature to be sig- 
nificantly suppressed with respect to what the T c of the 
hypothetically "clean" material would be. The dominat- 
ing in-plane disorder in as-grown crystals is likely to be of 
the same kind as introduced by the electron irradiation. 
Apart from unitary scatterers in the Cu02 planes, this 
also encompasses the "order parameter holes" induced 
by dopant oxygen and cation disorder in the rocksalt-likc 
layersj 18 ! 19 

The experimental data were confronted with a va- 
riety of theoretical models describing the reduction of 



c-axis superfluid density in terms of either quasiparti- 
cle dynamics or quantum fluctuations of the supercon- 
ducting order parameter phase in the CuC>2 layers. We 
find that the quantum phase fluctuation description 6,7 
yields excellent agreement as to the experimentally ob- 
served similar temperature dependences of A a & and A c , 
and quantitatively describes the temperature derivative 
d[fpi(T)/f P i(0)] 2 /d(T/T c ) 2 . However, it fails to de- 
scribe the dependence of the zero-temperature Josephson 
Plasma frequency on disorder strength. 

We therefore surmise that the reduction of f p i (0) with 
increasing disorder must be due to pair-breaking within 
the Cu02 layers. Data for X a b{T) and A C (T, T) are in 
agreement with scattering in the unitary limit by the 
irradiation-induced point defects in the Cu02-planes. 
Only one model consistently describes all aspects of the 
reduction of the c-axis superfluid density with tempera- 
ture and disorder strength. This is the Impurity Assisted 
Hopping model of Radtke et ai. j2I elaborated upon by Xi- 
ang and Wheatle y 23 ' 25 and by Kim and Carbotte»ii The 
model supposes a reduction of p c s {T) through hopping of 
nodal quasi-particles assisted by weak, highly anisotropic 
scattering by defects in the insulating SrO and BiO lay- 
ers. Candidates for such impurities are out-of-plane oxy- 
gen defects^ or cation disorder. From a fit of f p i (T) to 
the IAH model, we extract the energy scale A ~ 2.5fcsT c 
characterizing nodal quasiparticle excitations. This is re- 
markably close to the number obtained by Le Tacon et 
al. from anisotropic Raman scattering,— giving further 
confidence in the IAH interpretation. 
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